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The pasta phase in core-collapse supernova matter (finite temperatures and fixed proton fractions) 
is studied within relativistic mean field models. Three different calculations are used for comparison, 
the Thomas-Fermi (TF), the Coexisting Phases (CP) and the Compressible Liquid Drop (CLD) 
approximations. The effects of including light clusters in nuclear matter and the densities at which 
the transitions between pasta configurations and to uniform matter occur are also investigated. The 
free energy, pressure, entropy and chemical potentials in the range of particle number densities and 
temperatures expected to cover the pasta region are calculated. Finally, a comparison with a finite 
temperature Skyrme-Hartree-Fock calculation is drawn. 
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I. INTRODUCTION 

The complex structure of nuclear matter in the den¬ 
sity region approaching p s ~ 0.16 fm -3 (central density 
of heavy nuclei) at finite temperature (T < 20 MeV) crit¬ 
ically affects many astrophysical and nuclear physics phe¬ 
nomena. At low densities, the frustrated system called 
“pasta phase”, caused by the competition between the 
Coulomb interaction and the strong force, appears, and 
is constituted by different geometrical configurations, as 
the density increases M- 

The main interest in the pasta phase in core-collapse 
supernovae (CCSN) is that the neutrino opacity, which 
plays the main role in the development of a shock wave 
during the supernova collapse, is affected by its presence 
[2I. lid [Tlj|. At very low densities (up to 0.001 times the 
saturation density), light nuclei (deuterons, tritons, he- 
lions, a-particles) can appear |12l - tl7| . and, like the pasta 
phase, can modify the neutrino transport, which will 
have consequences in the cooling of the proto-neutron 
star Pilot. 

Following our previous works (see e.g. j9] and ref¬ 
erences therein), we use the Thomas-Fermi (TF) ap¬ 
proximation, where the surface effects are treated self- 
consistently, in the framework of relativistic mean field 
(RMF) models, to calculate the pasta phase for a fixed 
proton fraction, several temperatures and densities. Be¬ 
sides TF, we use the coexisting-phases (CP) method, 
which is numerically much faster than the TF one, and 
where the Gibbs equilibrium conditions are used to get 
the lowest free energy state, and the surface and Coulomb 
terms are added “by hand” (see e.g. 0)- Within this 
method, we also include the effect of light clusters. The 
compressible liquid drop (CLD) model is also considered. 
Unlike the CP approximation, this method takes into ac¬ 
count both the Coulomb and surface terms in the mini¬ 
mization of the total energy of the system. 

In the present work we perform a thermodynamical 
study of the phase transition between different shapes 
and we calculate the phase transformations occurring in 
the inner crust. We are interested in understanding the 
character of the crust-core transition and, for that reason, 


we study how the subsaturation instabilities are lifted by 
an appropriate description of the inner crust. In partic¬ 
ular, we want to identify the strong and weak points of 
each approach. 

We refer to phase transitions and phase transforma¬ 
tions as in Ref. [2(|, where the authors have stressed 
that there should be a distinction between “continuous 
transition” and “continuous transformation”. They make 
the difference between a transformation, which is a spe¬ 
cific path in the space of thermodynamic variables, and 
a phase transition, which is an anomaly of the thermo¬ 
dynamic potential considered in the total space of ther¬ 
modynamic variables. Having this in mind, an asym¬ 
metric system, containing a fixed proton fraction, is 
going through a first-order phase transition, since the 
first derivatives of the grand potential are discontinuous, 
though the thermodynamic transformations may result 
in a continuous evolution of the observables. 

We compare our results, obtained within the three 
specified methods, with the work of Raduta and Gul- 
minelli [21], Hempel and Schaffner-Bielich [22f , Zhang 
and Shen [23j], and Pais and Stone [|j. Raduta et 
al. developed a phenomenological statistical model for 
dilute star matter, in which free nucleons are treated 
within a mean-field approximation, and nuclei are con¬ 
sidered to form a loosely interacting cluster gas, with 
T = 1 — 20 MeV, y p = p p /p = 0 — 0.5 and p > 10 s 
g/cm 3 (p > 6.02 x 10~ 8 fm ~ 3 ), making it appropri¬ 
ate for CCSN description. They found that, for all sub¬ 
saturation densities, matter can be viewed as a contin¬ 
uous fluid mixture between free nucleons and massive 
nuclei. As a consequence, the equations of state (EoS) 
and the associated observables do not present any dis¬ 
continuity over the whole thermodynamic range and the 
expected first-order transition to uniform matter does 
not happen. Hempel et al. also used a statistical model 
for the EoS that can be applied to all densities relevant 
for supernova simulations, where an ensemble of nuclei 
and interacting nucleons are in statistical equilibrium. 
On the other hand, Zhang and Shen studied the non- 
uniform matter using a self-consistent TF approximation, 
though they only considered two geometrical configura¬ 
tions, droplets and bubbles. They have compared their 
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results with a parametrized TF calculation, where the 
surface energy and the nucleon distribution are calcu¬ 
lated differently, and they have reached the conclusion 
that the parametrized approximation is a reasonable one. 
We also compare our results with a 3D finite tempera¬ 
ture Skyrme-Hartree-Fock calculation [2d|, where four 
different Skyrrne interactions have been used, and where 
subtle variations in the low and high density transitions 
into and out of the pasta phase were found. 

The paper is organized as follows. In section El we 
briefly review the formalism used and in section IIIII the 
results are discussed. Finally, in section El some con¬ 
clusions are drawn. 


II. FORMALISM 

We consider a system of baryons, with mass M inter¬ 
acting with and through an isoscalar-scalar field (f> with 
mass m Sl an isoscalar-vector field V M with mass m v and 
an isovector-vector field with mass m p . When de¬ 
scribing npe matter we also include a system of electrons 
with mass m e . Protons and electrons interact through 
the electromagnetic field A A ‘. The Lagrangian density 
reads: 

c = y ] A + A + A + A + C p + A, 

i=p,n 

where the nucleon Lagrangian reads 

A = 4>i [ 7 /A-C^ - M*] A, (1) 

with 

iD* = id ** - g v V» - ^r • V 1 - (2) 

M* = M — g s (j) (3) 

and the electron Lagrangian is given by 

A = A bv (*<^ + eA^) -m e ]ip e . (4) 

The isoscalar part is associated with the scalar sigma 
(a) field </>, and the vector omega (w) field V p , whereas 
the isospin dependence comes from the isovector-vector 
rho (p) field b l p (where g stands for the four dimensional 
space-time indices and i the three-dimensional isospin 
direction index). The associated Lagrangians are: 

A = \ (d^d^ - m 2 s <p 2 - ^ A> 3 - 

A = + l -m 2 v VpV» + y g ^ v (VpVn 2 

C p = -±B„„ • B' ,u + l -m 2 p bp ■ 
r — _i p p^ v 

•‘-'7 — ^ l ll/ 


where = dpV v - d v Vp, B pu = dpb v - d v b p - g p (bp x 
bu ) and F pL , = 3 p A v Q v Ap. 

The model comprises the following parameters: three 
coupling constants g s , g v and g p of the mesons to the nu¬ 
cleons, the bare nucleon mass M, the electron mass m e , 
the masses of t he meso ns, the electromagnetic coupling 
constant e = -^/47 t/ 137 and the self-interacting coupling 
constants k, A and £. In this Lagrangian density, t is the 
isospin operator. 

We use the FSU parametrization |25|], expected to de¬ 
scribe well the crust [9|, even if it does not describe a 
2 M 0 neutron star. This parametrization also includes 
a nonlinear u>p coupling term, which affects the density 
dependence of the symmetry energy. This term is given 
by: 

C U] p = A v g 2 v g 2 bp-b^V fi Vr (5) 


The state that minimizes the energy of asymmet¬ 
ric nuclear matter is characterized by the distribution 
functions, /ofc±, of particles (+) and antiparticles (—) 
k = p,n, e, given by: 


/o ^= l +e te.)/T ' 3=P,n ^ 

with 

e 0j = \jp 2 + M * 2 , Vj = gj - g v V 0 {0) - ^ Tjb ( 0 0) (7) 

and 

/0e± = x + e(eoeTMe)/T , (8) 

with 

eoe = Vp 2 + ml, (9) 

where gk is the chemical potential of particle k = p, n , e. 

In the mean field approximation, the thermodynamic 
quantities of interest are given in terms of the meson 
fields, which are replaced by their constant expectation 
values. For homogeneous neutral nuclear matter, the en¬ 
ergy density, the entropy density, the free en ergy density, 
and the pressure are given, respectively, by [261 - 4281 ] : 
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£ = 


Y± J d PP 2e 0j {foj+ + hi-) 


m 


£>9 


v Vn+^V^ + 


v j/-4 , m P b 2 ■ ™ 2 


2 0 ' 2 

r2,2 

ptfv r 0 u 0 


^4>o + + 3As^V 0 2fe o> 
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+ 
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^ J dpp 2 [f 0j+ \nf 0j+ 

(! - foj+) ln(l - foj+) + foj- In hi- 
fi- - foj-) ln(l - foj-)\ > 

£-TS, 

j dp 7^ (/«+ + /«-) 




+ 2 *0 + 24 ^0 


-bi- 


- \<t >o - 7^0 + ^9 2 p9 2 vV£bl 


( 10 ) 


( 11 ) 

( 12 ) 


(13) 


Li is defined in Eq. © and, for all clusters, we have 
iD$ = id» - g vj V» e i±^, 

j = t,h,a,d. (21) 

The effective masses of the clusters are given by 

M* = 3 M - B t (22) 

Ml = 3 M - B h (23) 

M* =4 M - B a (24) 

M* d =2 M - B d (25) 

with their binding energies being B t = 8.482MeV, Bh = 
7.718MeV, B a = 28.296 MeV, and B d = 2.224 MeV. 
The fraction of a cluster, T), is given by Y) = pi/p,i = 
a, h,d,t. For the coupling constants, we consider g V j = 
Ajg v and g P j = \Zj— Nj\g p , where Aj is the mass number, 
and Z , Nj are the proton and neutron numbers, respec¬ 
tively [13 . The chemical potential of a cluster j is defined 
as pj = Njp n + Zjp p . More realistic parametrizations 
for the couplings of the light clusters have been proposed 
in [T 2 L flHj[. which should be implemented. 


For the electrons, we have 


£e 

7T 2 J 

f dpp 2 e 0e (/ 0e + + he-) , 


(14) 

Se 

1 

7 r 2 

J dpp 2 [foe+ In foe+ + (1 - 

- / 0e +) ln(l 

- he+) 


+ he- 

In/oe- + (1 - he-) ln(l - 

- he-)] ) 

(15) 


= £e~ 

TS e , 


(16) 

Pe 

1 

37T 2 

[dp— (foe+ + he-) ■ 

J Coe 


(17) 


We consider matter with fixed proton fraction that is 
neutrino free, and hence the neutrino pressure and energy 
density are zero [26.]. 


A. Light clusters 


We use the same prescription as in E3 to include light 
clusters (d = 2 H, t = 3 H, a = 4 He, h = 3 He) in the model. 
The Lagrangian density becomes 

L = 'y ] Li -V L a +L d -\-L e + £<j -\-Luj+L p + . (18) 

i=p,n,t,h 


The a particles and the deuterons are described as in 


L a = ±(iD^ a )*(iD^</, a ) - ±41 (M* a ) 2 4 a , 

C d = ±(iDl(f>d - iDd&d)*(iD d p,(t>dv - iD dv 4>dp) 

-\#z 


(19) 


B. Thomas-Fermi and Coexisting Phases 
appr oximat ions 


We use the Thomas-Fermi (TF) approximation to de¬ 
scribe the nonuniform npe matter inside the Wigner-Seitz 
unit cell, that is taken to be a sphere, a cylinder or a slab 
in three, two, and one dimensions 0[26]. In this approx¬ 
imation, npe matter is assumed locally homogeneous and 
at each point its density is determined by the correspond¬ 
ing local Fermi momenta. In 3D we consider spherical 
symmetry, in 2D we assume axial symmetry around the 
z axis, and in ID reflexion symmetry is imposed. In the 
TF approximation, fields are assumed to vary slowly so 
that baryons can be treated as moving in locally constant 
fields at each point [4j. In this approximation, the sur¬ 
face effects are treated self-consistently. Quantities such 
as the energy and entropy densities are averaged over the 
cells. The free energy and pressure are calculated from 
these two thermodynamical functions, using the usual 
expressions, see e.g. Ref. (27). 

In the Coexisting Phases (CP) method, matter is orga¬ 
nized into separated regions of higher and lower density, 
the higher ones being the pasta phases, and the lower 
ones a background nucleon gas. The interface between 
these regions is sharp. Finite size effects are taken into 
account bv a surface and a Coulomb terms in the energy 
density (l7j |. 

By minimizing the sum £ sur f + Seoul with respect to 
the size of the droplet/bubble, rod/tube, or slab one gets 

i 

surf = ‘2‘E'Couli (26) 


with 

Zcoul = (e 2 7T$) 1/3 [aD(p I p - pdp 1 )] 2/3 , (27) 


( 20 ) 
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where a = / for droplets, rods, slabs and a = 1 — / for 
tubes and bubbles, / is the volume fraction of phase /, 
a is the surface energy coefficient and $ is given by 


2 -D, 


1-2/D 


$ = 


D -2 1 

a —1 —In a 

D +2 


1 n _ i q 
D+2 


,D = 2 


(28) 


The Gibbs equilibrium conditions are imposed to get 
the lowest energy state, and, for a temperature T = T 1 = 
T 11 , are written as 

Mn = Pn, (29) 

A V'p 5 

pi = pi I , 


where I and // label the high- and low-density phases, 
respectively. When clusters are present, there are equi¬ 
librium conditions for them too |17j |. 

The total free energy density, and the total proton den¬ 
sity of the system are given by 

F = fF I + ( 1 - f)F u + F e + £ surf + Econu (30) 

P P = Pe = y v p = fp\ + (1 - f)pp, (31) 

where F l ,i = I, II, is the free energy density of the ho¬ 
mogeneous neutral nuclear matter, given by Eq. CE2D. F e 
is given by Eq. (fTTH) . and £ sur f and Seoul are the sur¬ 
face and Coulomb energies, given by Eqs. m and (E3, 
respectively. 


C. Compressible Liquid Drop model 

In the Compressible Liquid Drop model j~i~8l . l29l - [dH . the 
equilibrium conditions of the system are derived from the 
minimization of the total free energy j29], including the 
surface and Coulomb terms. This minimization is done 
with respect to four variables: the size of the geometric 
configuration, r^, which gives, just like in the CP case, 
Eq. (I2S1) . the baryonic density in the high-density phase, 
p 1 , the proton density in the high-density phase, p p , and 
the volume fraction, /. The equilibrium conditions be¬ 
come: 


Pn ~ Pn > 

pi = pi 1 - 


£surf 


m-fM-d 1 )’ 


P 1 = P 11 - £ 


surf 



( 32 ) 


P 1 / 


fO--f)(p P 



Note that there is an extra term in both the pro¬ 
ton chemical potential and in the mechanical equilibrium 
conditions, as compared to the ones obtained in the CP 
approximation, Eqs. (Eai. These terms arise from the 
inclusion of the surface and Coulomb terms in the min¬ 
imization of the total energy. The Coulomb repulsion 


induces an extra positive term while the surface tension 
reduces the cluster internal pressure. 

The total pressure, and the total proton chemical po¬ 
tential of the system are given by 

Ptot = ppPp + pnPn + PePe ~ Pj (33) 

Pp = fp*p + (1 - f)Pp, ( 34 ) 

where F is the total free energy density, given by Eq. ®, 
and / is the volume fraction of phase /. 


III. RESULTS 


In the present section, we discuss how the nuclear 
liquid-gas instability, occurring at subsaturation densi¬ 
ties for asymmetric nuclear matter, is partially lifted by 
an adequate description of the inner crust, allowing for 
the appearance of nonhomogeneous phases. In particu¬ 
lar, we will compare several physical quantities obtained 
within a TF calculation, a CP approach, supposing a 
zero thickness surface, and the CLD model, where finite 
size effects are included in a consistent way, with the cor¬ 
responding quantities for homogeneous matter. We will 
also discuss the effect of the inclusion of light clusters in 
the calculation. 

The free energy per particle, the pressure, the proton, 
neutron and baryonic chemical potentials, and the en¬ 
tropy per particle of the inner crust, obtained within the 
approaches referred above, are plotted in the following 
figures as a function of density or chemical potential, for 
the FSU interaction, two temperatures T = 4 MeV and 
8 MeV, and the proton fraction y p = p p /p = 0.3. The re¬ 
sults are shown for homogeneous matter (red), CP (blue) 
calculations with (dashed) and without (solid line) clus¬ 
ters, CLD (green solid line) and TF (points) calculations. 

For reference and to help the discussion, we show in 
Table H the symmetric nuclear matter properties for all 
the models we are using in this study to compare with 
our calculations with the FSU interaction: another RMF 
parametrization, TW [ 32 j, with density-dependent cou¬ 
plings, and four Skyrrne interactions, SkM* j33|, SLy4 
[HI], NRAPR [ 35 }, and SQMC700 [ 36 ij, chosen based on 
their overall performance in modelling a wide variety of 
nuclear matter properties ( 37 ). 

In Figure[Q we show the free energy per particle, F/A , 
as a function of the density. As expected, F/A is lowered 
when nonhomogeneous matter is present, making these 
states more stable. A second effect is the disappearance 
of the negative curvature that the EOS of homogeneous 
matter presents below saturation density. This effect is 
present in all the three methods considered. 

The light clusters are only present for very small densi¬ 
ties, and will start melting for p > 0.001 fm~ 3 . However, 
their presence lowers the free energy of the homogenous 
matter EOS and of the CP calculation, as can be seen in 
the inset panels. A TF calculation, inclu ding light clus¬ 
ters, should also be performed (see e.g. [38j). The CP 
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TABLE I. Symmetric nuclear matter properties at satura¬ 
tion density po (energy per particle B/A, incompressibility 
K, symmetry energy E sy m and symmetry energy slope L) for 
the FSU parametrization, and five other parameter sets for 
comparison. All the quantities are in MeV, except for po, 
given in fm~ 3 . 


Model 

po B/A 

K 

K 

J -'sym 

L 

FSU 

0.148 -16.3 

230 

32.6 

60.5 

TW 

0.15 -16.3 

240 

33 

55 

NRAPR 

0.16 -15.85 

226 

33 

60 

SQMC700 

0.17 -15.49 

222 

33 

59 

SkM* 

0.16 -15.77 

217 

30 

46 

SLy4 

0.16 -15.97 

230 

32 

46 


approach, which does not take into account in a consis¬ 
tent way the surface tension and Coulomb energy, over¬ 
estimates the effect of the clusterization, mainly at low 
densities. This problem is solved within the CLD ap¬ 
proach, which, taking into account finite size effects in 
the phase equilibrium conditions, gives results closer to 
the TF calculation. However, there is still some overes¬ 
timation of the free energy reduction with respect to the 
TF calculation, possibly due to the approximate descrip¬ 
tion of the surface energy. For the surface tension, we 
take a parametrization which is a function of the pro¬ 
ton fraction and the temperature, and was fitted to the 
results obtained from a relativistic Thomas-Fermi calcu¬ 
lation for a semi-infinite slab Close to the crust-core 
transition, all approaches, TF, CP and CLD, give simi¬ 
lar results, and predict a first-order phase transition to 
uniform matter. A first order phase transition at the 
crust-core transition has also been obtained within other 
approaches fzil. Eul]. 

When the temperature is increased to 8 MeV, simi¬ 
lar conclusions are drawn, the main differences being a 
decrease of the free energy and the density range of the 
pasta. A remark is in order: the first two points of the 
free energy for the TF calculation are above the homo¬ 
geneous matter value. This should be explained by the 
fact that, in the TF, we considered different values for 
the masses of the protons and neutrons, and by the fact 
that, for high temperatures, the precision in the calcu¬ 
lation begins to be very critical. Let us, however, point 
out that 8 MeV is already a quite high temperature, and 
thermal fluctuations of the clusters rod-like or slab-like 
will destroy the Wigner-Seitz structures according to [3f|. 

In Figure[2j the total pressure for the FSU interaction 
is plotted as a function of the baryonic chemical potential, 
ps, which is defined as in Ref. f22l |. 

PB = (1 Vp)Pn V VpiPp V (35) 

since we are performing a calculation with a fixed proton 
fraction and charge neutrality. The pressure is a smooth 
function of ps and, within the CLD and TF approaches, 
does not show any discontinuity. The CP approach does 
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FIG. 1. (Color online) Free energy per baryon as a function 
of the density for the FSU interaction, y v = 0.3, T = 4 MeV 
(top) and 8 MeV (bottom panel). Results with pasta (within 
CP, CLD and TF approaches) and for homogeneous matter, 
and including (for homogeneous matter and pasta within CP) 
or not clusters, are shown. The effect of these aggregates are 
only seen for very small densities (inset panels). 


not show any discontinuities between the intermediate 
shapes, however, at the onset of the pasta phase and at 
the crust-core transition, the pasta phase pressure does 
not match smoothly into the homogenous matter pres¬ 
sure. This is probably due to the non-consistent treat¬ 
ment of the surface energy. The free energy of the pasta 
matches the homogeneous matter free energy with a dif¬ 
ferent slope, both at low densities and at the crust-core 
density, giving rise to discontinuities in the pressure. The 
explicit inclusion of non-homogeneous matter, also true 
for the light clusters, increases the total pressure turning 
matter more stable. Consistently with the conclusions 
drawn for the total free energy, the contribution calcu¬ 
lated within the CP approach is larger than the contribu¬ 
tion evaluated within TF, showing the limitations of the 
first method. The CLD approach presents results very 
similar to the TF calculation. It is interesting to notice 
that the effect of light clusters is to increase the pressure, 
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FIG. 2. (Color online) Total pressure as a function of the 
baryonic chemical potential for the FSU interaction, y v = 
0.3, T = 4 MeV (top) and 8 MeV (bottom panel). Results 
with pasta (within CP, CLD and TF approaches) and for 
homogeneous matter, and including (for homogeneous matter 
and pasta within CP) or not clusters, are shown. 

making the system more stable. 

Let us now analyse the behavior of the pressure with 
the baryonic density. In particular, we will discuss how 
the instability at subsaturation densities in asymmetric 
nuclear matter is lifted in stellar matter. In Figure [3] 
the total pressure is plotted against the density, for the 
FSU interaction. First, it is interesting to notice that 
there are no instabilities, that is, there is no range of 
densities with a negative slope, except for small discon¬ 
tinuities, when the transition from one shape to another 
occurs, in the TF and CLD calculations. In the CP cal¬ 
culation, those discontinuities do not appear, because of 
the method itself: the Gibbs conditions, P 1 = P 11 , are 
imposed and no contribution from the surface tension is 
taken into account. There is, however, a jump at the on¬ 
set and melting of the inner crust, precisely due to the 
simplified description of the surface. The occurrence of 
discontinuities between the different pasta phases in the 
TF and CLD calculations originates from the simplified 
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FIG. 3. (Color online) Total pressure as a function of the 
density for the FSU interaction, y v = 0.3, T = 4 (top) and 8 
(bottom) MeV. Results with and without pasta, and including 
or not clusters, are shown. The effect of these aggregates are 
only seen for very small densities (inset panels). 


approach we have considered to describe the pasta phase, 
when only some geometries are included. However, in 
the TF calculation, even with this restriction, not always 
a transition between different shapes is discontinuous, e. 
g. the transition droplet-rod that occurs at pd- r = 0.024 
fm -3 , for T = 4 MeV, and at pd- r = 0.023, for T = 8 
MeV, is continuous. 

At the crust core transition we observe a discontinuity 
in all methods, with a larger density jump for the CP 
method and a smaller one for TF. A similar behavior has 
been obtained in other approaches @, [H, [22j, |4(| • Since 
the pressure does not change continuously into the ho¬ 
mogeneous matter, a Maxwell construction may be used 
to describe the transition from the non-uniform to uni¬ 
form matter, imposing a fixed proton fraction and charge 
neutrality, and equal pressures and baryonic chemical po¬ 
tential as defined in Eq. in both phases. Another 

aspect that should also be pointed out is that the onset 
of the pasta phases within the TF approach occurs at a 
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larger density than the CLD one. The TF approach has 
also several limitations, and a calculation with the ex¬ 
tended TF should be performed. However, we could also 
consider that, at smaller densities, the clustering onset is 
defined by the appearance of light clusters. 
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FIG. 4. (Color online) Fractions of the nucleons and clusters 
(top) and chemical potentials of the nucleons (bottom) as a 
function of the density for the FSU interaction, y v = 0.3 and 
T = 4 MeV, for homogeneous matter with clusters. 

We can also notice the effect of the clusters if we look 
at the inset panels in both plots. We observe that, for 
T = 4 MeV and p < 0.001 fm” 3 , the clusters slightly 
lower the pressure, bringing the homogeneous matter re¬ 
sult closer to the TF calculation. Above p > 0.001 fm” 3 , 
the clusters increase the pressure. Increasing the tem¬ 
perature, the clusters increase slightly the pressure in all 
range of densities shown, although for very low densities, 
the inclusion of light clusters still lowers the pressure. 
This occurs because the formation of light clusters in¬ 
creases the neutron fraction of the homogeneous matter, 
since light clusters are preferentially symmetric particles 
(deuterons or a particles) and, at high temperatures, the 
extra binding due to the cluster formation does not com¬ 
pensate the extra repulsion homogeneous matter with a 
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FIG. 5. (Color online) Pressure as a function of the density 
for the FSU interaction, y p = 0.3, and T = 4 MeV (top) and 
8 MeV (bottom panel). Homogeneous matter (HM) (solid), 
CP (dashed line), CLD (dash-point), and TF (points) results 
are shown for the baryonic (red), and electronic (blue) com¬ 
ponents. 

larger neutron fraction experiences. This is clearly seen 
in Fig. 21 where we plot, for T = 4 MeV, the neutron, 
proton, and light clusters fractions versus baryonic den¬ 
sity in the top panel, and the neutron and proton chem¬ 
ical potentials in the bottom panel. An increase of the 
neutron chemical potential above the homogeneous mat¬ 
ter one occurs close to the maximum of a clusters. This 
comes together with a decrease of the proton and neutron 
fractions, and a decrease of the proton fraction in the ho¬ 
mogeneous matter. Let us point out, however, that the 
thermodynamical function we should look at to discuss 
the extra stability the light clusters give to the system is 
the free energy, when the baryonic density is taken as a 
state variable. 

In order to understand the role of the electronic contri¬ 
bution to the total pressure, we plot in Fig. [5] the bary¬ 
onic (red) and electronic (blue) pressure components as a 
function of the density. This separation is not totally pos¬ 
sible for the TF calculation due to the coupling between 
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protons and electrons, induced by the Coulomb interac¬ 
tion. In the CLD method, and because of this separation 
problem, we are calculating P j,, given by Eq. m , as 
Pb = fP 1 + (1 — f)P 11 , unlike the other pressure plots, 
where it is calculated from the thermodynamical expres¬ 
sion, given by Eq. ®. The electronic pressure, P e , is 
given by Eq. m ■ Within the CP and CLD approaches, 
the electron contribution coincides with the homogeneous 
matter contribution. We see that for CP, the baryonic 
pressure is always positive because it corresponds to the 
pressure at the binodal surface, that is always positive. 
In the TF calculation, the baryonic pressure still exhibits 
a region with a negative incompressibility, unlike the to¬ 
tal pressure shown in Fig. [2] In fact, within TF, the 
negative slope of the pressure for nuclear homogeneous 
matter is only partially removed with the inclusion of 
droplets and the pasta phases. The electrons are respon¬ 
sible for making the total pressure positive, and with 
positive slope. It is also seen that the effect of clusteri¬ 
zation is washed out with the temperature, so that, for 
T = 8 MeV, the pressure for the pasta calculation does 
not differ much from the pressure for homogeneous mat¬ 
ter. Similar results were discussed in Ref. [2l| . where 
subsaturation densities are described within a statistical 
model, which considers a continuous fluid mixture of free 
nucleons and finite nuclei, the only difference being the 
absence of discontinuities between shape transitions. The 
discontinuity occurring in the electron pressure, within 
the TF calculation, is due to the normalization of the 
Coulomb field, which is set to zero at the Wigner-Seitz 
cell border. This gives rise to a jump going from slabs to 
tubes because, for the first geometry, the proton density 
is almost zero at r = Rws , while for the tubes it takes 
the largest value. 

In Figures [G] and 0 the neutron and proton chemical 
potentials are plotted as a function of the density, for 
the FSU interaction. Both quantities show an instability 
for homogeneous matter, explicitly seen through the pres¬ 
ence of a backbending: both chemical potentials decrease 
in a given range of densities, although the corresponding 
particle density is increasing. 

In Fig. [6l we plot the neutron chemical potential. The 
pasta phase practically removes the regions of instability. 
However, discontinuities are observed for all the interme¬ 
diate shape transitions in the TF and CLD calculations, 
and a large jump occurs at the onset and melting of the 
pasta phase in the CP approach, and only at the melt¬ 
ing of the crust, for the TF and CLD calculations. It is 
also seen that the light clusters are affecting the neutron 
chemical potential, particularly for T = 4 MeV: below 
p = O.OOlfm -3 , their presence decreases the chemical po¬ 
tential, while above that density, corresponding to the 
onset of the pasta phase, light clusters increase p n , as al¬ 
ready discussed in Fig. [3] This is due to the the appear¬ 
ance of more symmetric light clusters at larger densities, 
that induces a more neutron rich homogeneous matter. 

The proton chemical potential (Fig. [7]) behaves differ¬ 
ently: even though the back bending is reduced, it is still 



P (fm 3 ) 



FIG. 6. (Color online) Neutron chemical potential as a func¬ 
tion of the density for the FSU interaction, y v — 0.3, and 
T = 4 MeV (top) and 8 MeV (bottom panels). Results with 
and without pasta, and including or not clusters, are shown. 
The effect of these aggregates are only seen below 0.01 fm -3 . 

clearly seen. In homogeneous matter, for a given bary¬ 
onic density, the proton fraction is considered fixed and 
the charge neutrality is imposed. Consequently, there 
is only one global chemical potential associated to the 
global conservation of the baryonic number, Eq. (13511 . 
which is obtained from the derivative of the free energy 
with respect to the baryonic density. However, the cal¬ 
culation of the pasta phases in the present discussion 
was done imposing a global proton fraction together with 
the charge neutrality condition. Therefore, the neutron 
density and proton densities are independent, since neu¬ 
trons and protons can be exchanged freely between the 
two phases. Taking as degrees of freedom the baryon 
number and charge number, the neutron chemical po¬ 
tential is p n = pb and the proton chemical potential is 
p p = Pb+Pc- In neutrino free matter, the leptonic chem¬ 
ical potential is zero and the electron chemical potential 
is p e = — pc • Putting together the proton and electron 
chemical potentials, we conclude that p p + p e = Pb, see 
also the discussion in Ref. J4l|. In Fig. [5J the chemi- 
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FIG. 7. (Color online) Proton chemical potential as a function 
of the density for the FSU interaction, y v = 0.3, and T — 4 
MeV (top) and 8 MeV (bottom panels). Results with and 
without pasta, and including or not clusters, are shown. The 
effect of these aggregates are only seen below 0.01 fm -3 . 

cal potential /r p + fi e is plotted and no backbending is 
observed for the pasta calculations, showing that the in¬ 
stability has been totally raised, just like we have seen 
above for the neutron chemical potential. The large jump 
occurring at the slab-tube transition in Fig. [7] is related 
with the boundary conditions on the Coulomb field, when 
integrating the equation of motion. Measurable physical 
quantities are not affected by the boundary condition, 
but quantities, such as the chemical potential, are par¬ 
ticularly sensitive to the choice. 

In Figure [9j we plot the baryonic chemical potential 
/is, defined in Eq. (ESI, as a function of the baryonic 
density. For both temperatures shown, /is is a mono- 
tonically increasing function of the density, even for ho¬ 
mogeneous matter. However, by including a phase of 
non-homogeneous matter, the negative curvature of the 
chemical potential is removed. 

At the crust-core transition, TF, CP, and CLD give 
similar results. However, at the onset of the pasta phase, 
while TF and CLD link continuously to homogeneous 
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FIG. 8. (Color online) The /i p + /t e chemical potential as 
a function of the density for the FSU interaction, y p = 0.3, 
and T = 4 MeV (top) and 8 MeV (bottom). Results with 
and without pasta, and including or not clusters, are shown. 
The effect of these aggregates are only seen for very small 
densities. 


matter with clusters, the CP calculation presents a very 
large discontinuity, reflecting the non-consistent inclusion 
of the surface tension. This behavior is seen in all chemi¬ 
cal potential figures (except for fj, p at T = 8 MeV for the 
TF calculation). Light clusters reduce the homogeneous 
matter and the CP proton chemical potential, because 
clusters bring extra binding to the system, in particular, 
the a particles (see Fig. 0]). 

The discontinuities on the chemical potentials between 
the different geometrical configurations are an indication 
of the limitations of the present approaches. A rearrange¬ 
ment of matter and charge will wash out these disconti¬ 
nuities. Surface effects such as electrical double layers 
as the ones occurring on the boundary between charged 
solids and liquids, with an adsorption layer and a screen¬ 
ing layer, would give rise to a continuous behavior of the 
charge chemical potential. Also the Wigner Seitz approx¬ 
imation may disfavor an optimal matter rearrangement. 
Quantum molecular dynamics calculations go beyond the 
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FIG. 9. (Color online) Baryonic chemical potential as a func¬ 
tion of the density for the FSU interaction, y v = 0.3, and 
T = 4 MeV (top) and 8 MeV (bottom). Results with and 
without pasta, and including or not clusters, are shown. The 
effect of these aggregates are only seen for very small densi¬ 
ties. 


Wigner Seitz approximation and do not assume any spe¬ 
cific non-uniform structure of baryon matter, but con¬ 
sider an uniform background of electrons [5|-l7|. Using 
a large enough cell to include several units of the pasta 
structures, it has been shown recently in [4(| EU that a 
fee crystalline structure for the droplet phase would be 
favored with respect to the bcc one for densities just be¬ 
fore the transition to a rod like phase. An important 
conclusion of both approaches was the natural appear¬ 
ance of the typical pasta phases with rod, slab, tube, 
and bubble, in addition to spherical droplets, when no 
assumption on the structures was used. 

The behavior described above for both the proton 
and neutron chemical potential is also obtained in j23j . 
within a self-consistent TF description of the pasta phase, 
including only two configurations, droplets and bub¬ 
bles. The discontinuities obtained are attributed to the 
Coulomb field. Small discontinuities at low temperatures 
are obtained in [22j, within a statistical model, due to a 


stepwise increase of the average number of nucleons and 
protons of the clusters with the density. Although the 
overall behavior is the same in (2l| , where also a statisti¬ 
cal model is applied, and inhomogeneities are described 
as a continuous mixture of loosely interacting clusters, 
the chemical potentials show no discontinuities. The 
smooth increase of the average fragment mass fraction 
and cluster mass fraction explains this difference. 



P (fm -3 ) 
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FIG. 10. (Color online) Total entropy per baryon as a function 
of the density for the FSU interaction, y v = 0.3, and T — 4 
MeV (top) and 8 MeV (bottom panels). Results with and 
without pasta, and including or not clusters, are shown. The 
effect of these aggregates are only seen for very small densities. 

In Fig. [lOl the total entropy per baryon is plotted as a 
function of the density. The TF, CP and CLD calcula¬ 
tions lower the entropy per particle due to the formation 
of heavy clusters. For T = 4 MeV, and at low densi¬ 
ties, the same effect is seen in homogeneous matter with 
light clusters. On the other hand, in the CP calcula¬ 
tion, the inclusion of light clusters increases slightly the 
entropy because they reduce the formation of heavy clus¬ 
ters. When temperature is increased, the entropy for the 
pasta calculation gets closer to the homogeneous matter 
result because the nucleons drip out of the heavy clusters. 
The reduction of the entropy with clusterization was pre- 
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FIG. 11. (Color online) Crust-core transition densities, nor¬ 
malized to the nuclear saturation density po, as a function of 
the temperature, for the 3D-HFEOS (upward triangles) calcu¬ 
lation and the SkM* (green), SLy4 (pink), SQMC700 (blue), 
NRAPR (yellow) models, and the Thomas-Fermi (diamonds) 
calculation for the FSU (red) and TW (grey) models, and the 
CP (circles), and CLD (squares) results for the FSU (red) 
model. The proton fraction is fixed to 0.3. 


viously discussed in several works [HI, [2ll423j . Looking 
at Fig. 16 of Ref. [22} and Fig. 21 of Ref. [2lj], ap¬ 
plying a statistical description, a smooth decrease of the 
entropy per baryon with the density is obtained. Similar 
results have been obtained in [23} (see Fig. 2), within 
a self-consistent TF (STF) calculation for droplets and 
bubbles. 

In Fig. EH we show the crust-core transition densi¬ 
ties for the range of temperatures 4 to 10 MeV, and we 
compare them with the results found for the transition 
densities to homogeneous matter with a 3D finite temper¬ 
ature Skyrme-Hartree-Fock (3DHFEOS) calculation Q, 
and with the results found in Ref. j26j, for the TF calcula¬ 
tion with the density-dependent TW interaction. We see 
that the three approaches, TF, CLD and CP, give very 
similar results for the FSU interaction, and lower transi¬ 
tion densities than the values found within a 3DHFEOS 
calculation, done in the framework of Skyrme interac¬ 
tions. The TW parameter set has a higher transition 
density between shapes, and to uniform matter, than the 
FSU interaction, and very close to the values found for 
the Skyrme forces. This might be explained by the fact 
that the behavior of the symmetry energy for RMF mod¬ 
els with density-dependent couplings is much closer to 
the behavior of Skyrme forces [43J. It is known that the 
crust-core transition density decreases with an increasing 
slope L [44|, [4^] , and, therefore, this may influence the re¬ 
sults wtih SkM* and SLy4, having L = 46 MeV, while 
for all the other models, L ~ 60 MeV. Comparing models 
with a similar L, we conclude that, within the 3DHFEOS, 
the crust-core transition densities are larger by ~ 0.04po 
than the transition densities for FSU, and this difference 


decreases with increasing temperature. The 3DHFEOS 
calculation allows a larger freedom in the minimization 
of the free energy. However, differences may also be due 
to the different energy density functionals generated by 
each model. 
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FIG. 12. (Color online) Pasta phases for the FSU interaction 
within the TF (top), CP (middle) and CLD (bottom) calcu¬ 
lations for several temperatures. The proton fraction is fixed 
to 0.3. 
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TABLE II. Transition densities in the pasta phase for the TF, CP and CLD calculations for the FSU interaction, phm is the 
onset density of homogeneous matter. These values are the ones represented in Fig. 1121 
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TF CP CLD 

TF CP CLD 
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0.024 
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0.039 

0.052 0.085 - 
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6 

0.025 
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0.040 
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0.068 - 

0.087 0.086 

0.088 

8 

0.027 

0.032 

0.025 

0.040 0.048 

0.040 

- 

0.061 - 

0.079 0.078 

0.079 
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0.031 

0.032 

- 

- 0.047 

- 

- 

0.047 - 

0.068 0.066 
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In Ref. [23j, even though only droplets and bubbles 
were considered for simplicity, the transition densities to 
homogeneous matter are very close to our results with 
T = 10 MeV (see e.g. Fig. 1 of Ref. @). 

Let us now compare how the different approaches de¬ 
scribe the transition between the different shapes. In Ta¬ 
ble EH we give the densities at the droplet-rod, rod-slab, 
slab-tube, tube-bubble transitions, and the onset density 
for homogeneous matter. 

These transition densities have been plotted in Fig.fl2l 
where the density range for each geometric configuration 
is shown for the TF (top panel), CP (middle panel) and 
CLD (bottom panel) calculations, and a range of temper¬ 
atures 4 to 10 MeV. It is interesting to notice that in the 
CLD approach, tubes and bubbles do not exist, and at 
T = 10 MeV, the pasta geometries no longer exist. For 
the CP calculation, the tubes and bubbles are also not 
favored: droplets, rods and slabs are present for all tem¬ 
peratures, but tubes only briefly appear, at T = 4 and 6 
MeV, and no bubble configuration was found within this 
calculation. For the TF approximation, we found all the 
shapes, except for the slabs, that only appear at T = 4 
MeV. 
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FIG. 13. Pasta phases for the TW interaction within the TF 
calculation for several temperatures. The proton fraction is 
fixed to 0.3. (data taken from El). 


However, it is important to stress that the appearance 


and density range of the different geometries is model 
dependent, and sensitive to properties such as the sym¬ 
metry energy. In particular, in Figure [T3l where the dif¬ 
ferent pasta geometries are represented for the density 
dependent RMF model TW [32[, for a proton fraction of 
0.3, we can see that all the five geometrical configurations 
are present until T = 8 MeV. The transition densities to 
uniform matter, as we already saw, are represented in 
Fig. [TTJ la Ref. [43], it has been shown that the be¬ 
havior of the symmetry energy, and its derivatives with 
the density, for relativistic nuclear models with density 
dependent couplings, is much closer to the behavior of 
Skyrme forces. 

In Ref. [||, it was discussed the influence of a correct 
treatment of the Coulomb interaction on the extension 
of each pasta geometry. In particular, in a calculation 
excluding the Coulomb field, and including the Coulomb 
energy, as was done in the CP calculation, the bubble 
geometry was not present with Y p = 0.1, and the slab 
configuration was found in a wider density range with 
Y p = 0.3. The same authors have also discussed the role 
of the surface tension, and have shown that a smaller sur¬ 
face tension favors a larger variety of geometries. In our 
CP and CLD calculations, the surface tensions have been 
calculated within the model, however the treatment is not 
selfconsistent, and for the larger densities, when the back¬ 
ground nucleon gas becomes denser, the surface energy 
used is probably too high and the neutron skin should 
have been included explicitly in the surface energy, as in 
Ref. [4f|. This also explains why the crust-core transi¬ 
tion density calculated within the CP approach is smaller 
than the TF result. 


IV. CONCLUSIONS 

In this work, we have studied the pasta geometries that 
appear in core-collapse supernova events, and in the in¬ 
ner crust of neutron stars, within three different approx¬ 
imations, the Thomas-Fermi, the Coexisting Phases and 
the Compressible Liquid Drop calculations, all within the 
single nucleus approximation. While the first is a selfcon¬ 
sistent calculation, where the Coulomb interaction and 
surface energy are adequately described, the other two 
approaches are non-selfconsistent. They use as an input 
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the surface energy for semi-infinite matter as a function 
of the proton fraction and temperature, obtained within 
a Thomas Fermi calculation. In the CLD approach, the 
equilibrium conditions between the liquid and gas phases 
take into account the Coulomb energy and surface con¬ 
tributions, contrary to the CP method. An improvement 
of the last two methods is tightly related with a realistic 
description of the surface of the clusters. 

We have introduced light clusters into our system to 
understand their effect on the EoS. We observed that 
their effect is only noticeable at very low densities, be¬ 
fore melting. It was shown that taking light clusters into 
account always lowers the free energy. The inclusion of 
light clusters in all the methods used in the present work 
to describe the inner-crust will allow going beyond the 
single nucleus approximation. 

We were also interested in characterizing the transi¬ 
tion to uniform matter. For that effect, we have plotted 
the free energy, pressure, entropy and chemical potentials 
to observe if there were any discontinuities. We realized 
that the density range of the pasta phase and the crust- 
core transition density decrease with increasing temper¬ 
ature, as expected, however the melting temperature 
of the different pasta phase geometries depends on the 
model properties. The stable geometries that we found 


depend on the parametrizations used, and the proper¬ 
ties that influence them should be investigated. Also, 
the jumps in the pressure and chemical potentials, as a 
function of the density, could indicate a first order phase 
transition to uniform matter. Within the CP method, the 
description of a non-homogeneous phase gives unrealis¬ 
tic results at low densities, though it predicts concordant 
transition densities to uniform matter. All the methods 
considered in this study show a very good agreement with 
respect to the transition density to homogeneous matter, 
and, in particular, the TF and CLD calculations give 
very similar results in the whole range of densities and 
temperatures considered. 
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